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Abstract. We present an analysis of neutrinos de- 
tected with the Baikal neutrino telescope NT200 for 
correlations with gamma-ray bursts (GRB). No neu- 
trino events correlated with GRB were observed. As- 
suming a Waxman-Bahcall spectrum, a neutrino flux 



upper limit of < 1.1 x 10 



'GeV 



was obtained. We also present the Green's Function 
fluence limit for this search, which extends two orders 
of magnitude beyond the energy range of the Super- 
Kamiokande limit. 

Keywords: Neutrino telescope, BAIKAL, Gamma- 
ray burst 

I. Introduction 

The Baikal neutrino telescope NT200 |[l], lIU is op- 
erating in Lake Baikal, Siberia, at a depth 1 . 1 km since 
April, 1998. NT200 consists of 8 strings of 70 m length: 
7 peripheral strings and a central one. Interstring dis- 
tances are about 20 m. Each string includes 24 pairwise 
arranged optical modules (OM). Each OM contains a 
37-cm diameter hybrid photodetector QUASAR-370. 

A number of relevant physics results has been ob- 
tained so far with the NT200 telescope, e.g. limits on 
the diffuse flux of extraterrestrial high energy neutrinos, 
limits on neutrino fluxes from Dark Matter annihilation 
(Sun, Earth), and on the flux of relativistic and slow 
magnetic monopoles (|2l, O, ||4 |. 

This work is devoted to the search of neutrino events 
correlated with observations of more than 300 gamma- 
ray bursts (GRBs) reported from 1998 to 2000 by the 
Burst and Transient Source Experiment (BATSE) [5 |. 



The detection strategy for neutrino events with the 
NT200 telescope is based on a search for Cherenkov 
light from relativistic up-going muons produced by 
neutrino interactions. Information about the GRB time 
and location on the sky allows to reduce the atmospheric 
muon background and, as a result, significantly increases 
the sensitivity of the neutrino telescope to neutrino 
events correlated with GRB. 



II. Experimental data 

For the present analysis, the experimental data ob- 
tained with NT200 from April 1998 to May 2000 were 
used. The selected data sample contains those events 
which were formally reconstructed as up-going muons. 
Taking into account the high level of background for 
directions close to the horizon, only events with zenith 
angles larger than 100° were selected. The average rate 
of such events was 0.037 Hz. Most fake events are due 
to misreconstructed muons close to horizon and to muon 
bundles. 

For the present analysis of time and directional corre- 
lations with NT200 events we used information about 
GRB location, time, duration Tgo, and location error 
from the basic BATSE 4B catalog (5] (triggered bursts) 
and from the catalog of non-triggered GRB |i6J. The 
error distribution of BATSE GRB locations was taken 
from |7|. A total of 303 GRBs (155 tilggered and 148 
non-triggered) at zenith angles larger than 100° and 
occurring during periods of stable operation of NT200 
have been selected. 
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III. Data selection criteria and detector 

EFFECTIVE AREA 

The optimization of the data selection criteria was per- 
formed on the basis of simulated neutrino events [8^] and 
events of atmospheric muon background ||9l in NT200. 
Taking into account the varying NT200 configurations 
during the considered time period, calculations have 
been performed for nine basic detector configurations, 
most of them closely corresponding to the real status of 
the detector 

The results of the reconstruction of simulated events 
were used to estimate the reconstruction efficiency and 
to calculate the background. Event reconstruction and 
data selection for NT200 are described in detail in |8|. 
There, selection criteria were designed and optimized 
for atmospheric neutrino separation. They provide a re- 
jection factor of atmospheric muons larger than 10^. For 
GRB, the aditional information about detection time and 
location on sky, however, allows softening the require- 
ments to the background rejection. This increases the 
registration probability for useful events and therefore 
greatly increases the effective neutrino detection area. 
Following the approach of [!8|, Phu x Pnohit and Zdist 
were chosen as basic parameters for event selection. 
Zdist is the maximal distance between all projections 
of the triggered OM coordinates onto the reconstructed 
muon trajectory. Phu is the normalized probability of 
fired channels to be hit, and Pnohit is the probability of 
non-fired channel not to be hit. 

For the present correlation analysis, two sets of criteria 
for event selection were chosen: 
Cut-A: 

{Zd^st > 30m)k{PMt X PnolM > 0.1)&(* < 10°), 

Cut-B: 

iZd^st > 30m)k{^ < 5°), 

were is the angle between up-going muon and GRB- 
direction. 

Cut-A dominantly selects neutrinos with energies be- 
low ~ 10^ GeV. Cut-B allows a significant extension 
of the energy range, but the expected background is 
approximately four times larger than for Cut-A. 

Calculating the effective area of NT200 for the two 
sets of criteria, we took into account the absorption 
of neutrinos passing through the Earth, as well as the 
production, propagation, detection and reconstruction 
within a given angular cut 5* of muons. The calculated 
effective areas for the Cut-A and Cut-B samples are 
presented in Fig. 1 as a function of neutrino energy. 

The effective areas for the two sets are close to 
each other up to 10^ GeV. For larger energies, the 
effective area for Cut-A stays essentially constant. The 
behavior for E > 10^ GeV is largely defined by neutrino 
absorption in the Earth. 

The energy range of the NT200 sensitivity was esti- 
mated for an E^"^ neutrino spectrum. The 90% sensitiv- 
ity range of NT200 extends up to - 10*^ GeV and ~ lO'^ 
GeV for selection criterion A and B, respectively. 
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Fig. 1. NT20() effective area averaged over zenith angles between 
100° and 180° as a function of neutrino energy, for selection criteria 
Cut-A and Cut-B. 

A global estimation of the background, expected for 
the GRB search time window, was obtained from the 
NT200 raw event rate (0.037 Hz), the calculated atmo- 
spheric muon rejection factor for Cut-A and Cut-B, and 
the total GRB duration Tgrb 1-8 x 10^ s). Tgrb 
was calculated as the sum of the Tgo intervals for all 
303 triggered and non-triggered GRB. For compensation 
of possible event time uncertainties, five seconds were 
added at both sides of the time intervals Tgo. For 
cases of missing information on Tgo (about 25% of 
the triggered GRB), a fixed time interval was used. 
The expected background values, obtained from the full 
detector livetime interval, are 0.8 and 3.1 events for Cut- 
A and Cut-B, respectively. 

IV. Data analysis and results 

The basic objective of this study was to verify a suf- 
ficiently precise detector responce simulation, to search 
for events correlated with GRB, and to provide a solid 
background estimation. 

To check the simulation procedures, the calculated 
atmospheric muon rejection factors were compared to 
experimental values. The results are presented in Table 
I for different criteria Phit x Pnohit (Zdist is set to 
30 m). The simulated values are in agreement with the 
experimental results within the systematic error of our 
calculation, about 20%. 

TABLE I 

Atmospheric Muon Rejection Factor: Experiment and 
Simulation 



Phit ^ Pnohit 


> 0.1 


> 0.2 


> 0.3 


Experiment 
Model 


0.053 
0.062 


0.012 
0.014 


0.0035 
0.0040 



This search for correlation with NT200 neutrino 
events uses 303 GRBs (triggered and non-triggered) 
selected from the total sample of 736 BATSE GRBs 
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recorded in 1998-2000. The search for events selected 
according to criteria A and B was done for the time 
intervals Tgrb (see section III). For a detailed back- 
ground estimation, we used a time interval ±1000 s with 
respect to the start time of the GRB (excluding the Tq rb 
window) and an angle between GRB and up-going muon 
* < 10°. 

No events were found according to criterion A and 
one event by criterion B. Table II shows the number of 
signal and background events, the event upper limit /igp 
that was obtained in accordance with [10], the number 
of GRB corrected to coefficient /3 (the probability to 
detect events associated with gamma-ray burst in the 
given angle interval given the location error of the 
GRB), and the 90% C.L. upper limit on the number of 
events per GRB A'go — /igo / {Ngrb x P)- Data are 
presented for all GRB (triggered and non-triggered), as 
well as separatly for triggered GRB selected according 
to Cut-B. 

TABLE II 
Results of GRB Analysis 



Selection 


Signal 


Backgr 


M90 


Ngrb x P 




Ciit-A, all 





0.56 


1.9 


236 


0.0085 


Cut-B. all 


1 


2.7 


2.1 


199 


0.010 


Cut-B, trig 


1 


1.6 


2.8 


120 


0.023 



No excess of events associated with a GRB was 
observed. The limit on the neutrino flux associated with 
gamma-ray bursts was obtained using the approach from 
ifTTll . According to this approach, the limit F{E^) is 
presented as a function on neutrino energy, the "Green's 
function": 



F{E,)^N^o/Seff{E,), 



(1) 



Seff{E^) is the detector effective area, A^go the 90% 
C.L. upper limit on number of events per GRB. The 
main advantage of this approach is that its result does 
not depend on assumptions about the neutrino energy 
spectrum. 

Figure 2 shows the 90% C.L. upper limits of the GRB 
neutrino fluence Green's function F{E^) for NT200 
(for Cut-B, all), Super-Kamiokande [11 J and AMANDA 
lfT2ll . The Super-Kamiokande and NT200 limits are 
mainly for GRB from the southern sky, while the 
AMANDA limit is for the northern sky. 

Since predictions for the energy spectrum of neutrinos 
from GRB differ from model to model, we prefer 
to present our basic experimental result as a Green's 
function F{E^), which allows to calculate limits for any 
neutrino energy spectrum. In addition, we translate our 
result to a benchmark spectrum. We have chosen that of 
E.Waxman and J.Bahcall [MJ, and E.Waxman |[l5l, lfT6l . 
Following these works, the muon neutrino differential 
flux ^w-b{Eu) in the energy range up to 10 PeV is: 

El <i>^-^(i?.) = A^-^ X min(l, E./E,,), (2) 
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Fig. 2. 90% C.L. upper limits on the GRB neutrino fluence Green's 
function F{E^) for NT200, Super-Kamiokande and AMANDA. 



-B^6=100 TeV, Aw-^B ~8x 

The Model Rejection Factor MRF for the Waxman- 
Bahcall spectrum was calculated as: 



MRF - Ngo/Ne 



(3) 



where A^go is the upper limit on the number of events per 
GRB and N^x the expected number of events, calculated 
for the given spectrum as: 



Earth 



[E,) X S,,ff{E,) X {ATi/n)dE,. (4) 



Here, n w 2.2 x 10 ^ is the average GRB rate in 
AiTsr (~ 700 events within the detection range of the 
BATSE per yeai-) and $f '^''*'*(i;^)=0.5 x ^"^-^{E^) 
the neutrino flux at the Earth. 

Taking into account that the estimation of the expected 
event number in the given approach is made for for the 
BATSE burst detection rate, the MRF was calculated 
only for triggered GRB {Ngrb x /3 = 120, see Table 
II). From that Green's function (approximately twice that 
of the (Cut-B, all) sample, see Fig. 2), the resulting MRF 
value is 2.8 x 10^, and the corresponding GRB neutrino 
flux limit is 



Ft < 1.1 X IQ-'^GeVcm' 



(5) 



This diffuse limit is considerably weaker that of the 
AMANDA muon analysis llT2l . and two times higher 
that their cascade analysis ifTSl . In view of a search for 
bright individual GRBs, our result may be considered 
as complementary to AMANDA since variations in ab- 
solute energy output, Lorentz factor and distance might 
lead to a GRB neutrino detection with a less sensitive 
detector, while that source was outside the other detector 
field of view. NT200H- is presently complementing the 
ANTARES detector, in particular for those cases where 
the GRB is above horizon with respect to this instrument. 
We also not that, normalized to a single GRB, NT200 
exceeds the sensitivity of Super-Kamiokande by a factor 
of 2 for neutrino energy above 1 TeV. 
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V. Conclusion 

We have presented results of a search for neutrino 
induced muons detected with Baikal Telescope NT200 in 
coincidence with 303 gamma-ray bursts recorded from 
1998 to 2000 by BATSE. NT200's field of view covers 
most part of the Southern hemisphere. No evidence 
for neutrino-induced muons from gamma-ray bursts is 
found. The resulting Green's Function fluence limit for 
this search extends that of Super-Kamiokande by two 
orders of magnitude in energy. Assuming a Waxman- 
Bahcall spectrum, a neutrino flux upper limit of 
El^v < 1-1 X lO^'^GeVmi^'^s^^sr^^ is obtained. 
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